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ABSTRACT 

Several recent studies have shown that the star cluster initial mass function (CIMF) can be 
well approximated by a power law, with indications for a steepening or truncation at high 
masses. This contribution considers the evolution of such a mass function due to cluster dis- 
ruption, with emphasis on the part of the mass function that is observable in the first ~1 Gyr. 
A Schechter type function is used for the CIMF, with a power law index of —2 at low masses 
and an exponential truncation at M* . Cluster disruption due to the tidal field of the host galaxy 
and encounters with giant molecular clouds flattens the low-mass end of the mass function, 
but there is always a part of the 'evolved Schechter function' that can be approximated by a 
power law with index —2. The mass range for which this holds depends on age, r, and shifts to 
higher masses roughly as t 0,6 . Mean cluster masses derived from luminosity limited samples 
increase with age very similarly due to the evolutionary fading of clusters. Empirical mass 
functions are, therefore, approximately power laws with index —2, or slightly steeper, at all 
ages. The results are illustrated by an application to the star cluster population of the interact- 
ing galaxy M51, which can be well described by a model with M* = (1.9 ± 0.5) x 10 5 M 
and a short (mass-dependent) disruption time destroying M„ clusters in roughly a Gyr. 

Key words: galaxies: star clusters - open clusters and assocations: general - globular clus- 
ters: general 



1 INTRODUCTION 

The survival chances of young (< 1 Gyr) star clusters and 
■ how this may, or may not, depend on their mass, M, and 
environment has been the topic of quite some debate re- 
cently. Two models for the early evolution of clusters have re- 
cently become available, with a very different role of M in 
the disruption description: a mass-depende nt disruption model 
( Boutlouko s"& Lamersl2003l hereafter BL 03: lLamers et all2005al : 
lLamers. Gie les & Portegies Zwart 200 5b|) and a mass-indep endent 



disruption model ( frail et al.l 120051 : IWhitmore et alj 120071) . Both 
disruption models assume a continuous power law cluster ini- 
tial mass function (CIMF) with index —2, based on empirical 
derivations of the mass function of young (~ 10 Myr) clusters with 
masses between ~1O 2 M0 and ~1 5 Mn in a v a riety o f galac- 
tic environments (see for example IZhang & Falll dl999|) for the 
case of the cluste rs in the Antennae galaxies; iBik et alj ( 120031) for 
clusters in M51; [McCradv & Graham (200% for cluster in M82; 



Ide Griis & Goodwirj (2008) for the SMC an d lde Griis etal] f2003a) 
for a compilation of clusters in different galaxies j]]. 

Based on a study of the age and mass distributions of star 
clusters in four different galaxies, BL03 conclude that the lifetime 
of star clusters, or disruption time-scale, t^ s , depends on M as 
tdis oc A/ 7 , with 7 ~ 0.62 ± 0.06 and with the proportionality 
constant dependent on galactic environment. This mass-dependent 
disruption (MDD) model is supported by our theoretical under- 
standing of the disruption of star clusters, since the same scaling 
of tdis with M, i.e. the same value of 7, was found from iV-body 
simulations of clusters dissolving und er the combined effect of in- 
ternal relaxation and exter nal tides {Baumgardt & Makino [20031 : 
iGieles & Ba umgardt 2008), and also for the disruption of clusters 
due to external perturbations, or 'shocks'. In fact, tdis due to shocks 
depends on both Al and the half-ma ss radius, n,, s i nce tdis scales 
with the cluster density within n, (Spitze il l 19581 : lOstriker et all 
fl972h . which combined with the shallow dependence of n, on M 
for (young) clusters (rh oc M 01 .e.g. lLarsenl2004l) . leads to a sim- 
ilar value of 7 as for tidal disruption. 



1 The observationally derived CIMF from young (<10 Myr) clusters will 
be referred to as the CIMF cmp . This is because the number of massive clus- 
*E-mail: mgieles@eso.org ters younger than ~ 10 Myr is not necessarily sufficient to reveal the shape 

of the CIMF at high masses. See Section |2"TTl for more details. 
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MDD flattens the low-mass end of the cluster mass function, 
dN/dhl, at old ages. iFall & Zhand j200 lh show that for a con- 
stant mass-loss rate (7 = 1), the mass function of a single-age 
cluster population evolves towards a flat distribution at the low- 
mass end, independent of the initial shape of the mass function. 
Such a mass function is in good agreement with that of the Galactic 
globular cluster system. However, the globular cluster mass func- 
tion (GCMF) is usually derived from th e luminosity funct i on us- 
ing a constant mass-to-light ratio, M/L. Man dushev et alj J 1 99 lh 
find from dynamical mass estimates of Galactic globular clusters 
that their M/L is actually an increasing function of M, which 
could be the result of the p referential depletion of low-mass stars 
(Kruii ssen & Lamersll2008l) . Such a mass-dependent M/L alters 
the shape of the GCMF and could allow for smaller values of 7 
(Kruijssen & Portegies Zwart, in prep). A general expre ssion for the 
logarit hmic slope at low masses due to MDD is 7 — 1 dLamers et al.l 
l2005ah . 

An MDD model for the early evolution of star clusters was 
introduced by BL03. In the four galaxies studied by BL03, only 
masses for clusters in M33 and M5 1 were available. The mass func- 
tions of clusters in these galaxies were presented for all ages, so 
no evidence for a flattening of the mass function for a selected set 
of old clusters was provided, although these integrated mass func- 
tion s are slightly flatter at low masses than a power law with index 
— 2. lde Griis et alj d2003bh report a turnover in the mass function of 
clusters in region B of M82 and they claim that the population is 
approximately coeval (~1 Gyr). This result was later refuted, since 
from spectr oscopic studies a larger age r ange with a younger mean 
was found dKonstantopoulos et al.ll2008l) and the turnover was at- 
tributed to detection incompleteness d u e to the extent of clusters 
dSmith et alj|2007l : iMavva et alj|2008h - Ide Griis & Anders! d2006l) 
determine the mass function of LMC clusters in different age bins 
and find that it steepens with age, from a power law function with 
an inde x of —1.8 at ^ 1 Myr to one with an index of —2.2 at 
~1 Gyr. IPiskunov et all d2008t) also report a steeper mass func- 
tion at older ages for Milky Way open clusters. This steepening 
is contrary to what is expected from the evolutio n of a continuous 
power law C IMF with mass-dependent disruption. Goudfro oii et all 
d2004l2007t) show that the intermediate age (few Gyrs) cluster pop- 
ulations of the merger remnants NGC 1316 and NGC 3610, respec- 
tively, provide evidence for a turnover in their mass function, based 
on a flattening of the cluster luminosity function at low luminosi- 
ties. But these clusters have already evolved over a significant frac- 
tion of a Hubble time, so, it is not at all convincing that the mass 
function of clusters younger than a Gyr gets shallower than a power 
law with index —2. 

IZhang & F all (1999) determined the mass function of clus- 
ter in the Antennae aged between 25 and 160 Myr and show that 
it has essentially the same shape as that for the cluster s between 
2.5 an d 6 Myr, namely a power law with index —2. IFall et all 
(2005) show that the age distribution, dN/dr, of mass limited sub- 
samples of clusters in the Antennae galaxies declines as r -1 . Such 
an age distribution can result from a constant formation history 
combined with a 90% disruption fraction each age dex. Together 
with the similarity of the mass functions at different ages the au- 
thors concluded that the disruption of star clusters in the first Gyr 
is independent of their mass. Several claims have been made that 
this mass-independent disruption (MID) model also describes the 
age distribution of clusters in other galaxies, such as the SMC 



dChandar et alll2006l) and M33 dSaraiedini & ManconelkoOTh . but 
these results have later been ascribed to detection incompleteness, 
since luminosity limited cluster samples, not affected by disrup- 
tion, also have an age distribution that sca les approximately as r _1 
dGieles. Lamers & Portegie s Zwartl2007l) . 

The consequence of the 90% MID model is that the number of 
clusters in logarithmic age bins is constant, which for a continuous 
power la w CIMF results in a cons tant maximum cluster mass in 
such bins. Gieles & Bastianj d2008h showed that this is indeed the 
case for the Antennae, but through a comparison to cluster popu- 
lations in six other galaxies, they also show that the cluster popu- 
lation of the Antennae galaxies is unique in that sense. The MID 
model fails to reproduce the first few 100 Myr of the (mass limited) 
dN /dr in other galaxies ( see the discussion in lGieles et alj|2007l 
an d iGieles & Bastianll2008l) , so it seems not to be a 'universal' sce- 
nario for cluster evolution, as claimed by Whitm ore et alj {2007). 

In summary, the MDD model can not reproduce the correct 
mass functions at old ages and the MID model can not reproduce 
the correct age distribution at young ages in galaxies other than the 
Antennae. 

This study searches for the explanation of this problem by 
abandoning an assumption that both models make, na mely a con- 
tinuou s power law CIMF. Alternatively, a truncated ISchechteil 
(1976) function (equation [TJ is considered for the parent distri- 
bution from which cluster masses are drawn and this function is 
evolved with mass-dependent disruption. In Section [2] arguments 
for the choice of this truncated distribution function are given and 
its basic properties are presented. In Section [3] an analytical model 
for evolved Schechter mass functions is presented and a compar- 
ison to the cluster population of M51 is provided in Section [4] A 
discussion and conclusions are given in Section|5] 



2 A SCHECHTER FUNCTION FOR THE CLUSTER 
INITIAL MASS FUNCTION 

2.1 The need for a truncation 

Assum e that initial cluster masses, Mi, are drawn from a lSchechterl 
dl976l) parent distribution function of the form 

^=AM- 2 exp(-M i /M»), (1) 

where Af* is the mass where the exponential drop occurs and A 
is a constant that scales with the cluster formation rate (CFR). The 
constant A can also be taken as a function of time, such that equa- 
tion Q} describes the cluster formation history and thus represents 
the probability of a cluster with an initial mass between Mi and 
Mi + dMi forming at a time between t and t + dt. Since A will 
be assumed to be constant in most cases throughout this study and 
since the focus will be on the evolution of the mass function the no- 
tation dN/dMi is preferred, though dN/(dA4idt) would be more 
precise. 

As mentioned in footnote 1, this distribution function for the 
cluster initial mass function (CIMF) is not necessarily reflected 
from the empirically derived CIMF (CIMF cmp ). This is because 
usually only clusters younger than ~ 10 Myr are used to determine 
the CIMFcmp. This time interval is short compared to the range of 
cluster ages in a typical cluster population (few Gyrs). When sam- 
pling only a low number of clusters from the distribution function 
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in equation |QJ, the most massive cluster actually formed, M max , 
can be less massive than M*, and then the truncation is not recog- 
nised from the CIMF cmp . The analytical form of equation Q} was 
proposed for the luminosity function galaxie s, where L« (inst ead of 
M*) is the characteristic galaxy luminosity dSchechteill 19761) . This 
functional form follows from a stochastic self-si milar model for the 
origin of galaxies from self-gravitating gas (e.g. lPress & SchechteJ 
1 19741) . On stellar scales a power law (self-similar) mass spectrum at 
low masses with an exponential cut-off at high masses also follows 
from t heoretical models of fr agmentation in a collapsing molecular 
cloud dSilk&Takahashill979h . For a constant cluster formation ef- 
ficiency the CIMF reflects the shape of the giant molecular clouds 
(GMCs) and cores from which the clusters form. Collisions be- 
tween cores within the GMCs and destruction of the most massive 
cores by the star formation process also le ads to a Schechter type 
(equi librium) mass function for clusters fMcLaughli n & Pudritj 
[1991) . 

The choice for a truncation is also motivated by several indica- 
tions that the high-mass end of the cluster mass function is steeper 
when a larger age range, and thus a larger total number of clusters, 
is considered. The most important arguments are discussed below. 



2.1.1 The most massive cluster in logarithmic age bins 

iHunter et alj d2003l) studied the evolution of the most massive 
cluster mass, M max , in equally spaced logarithmic age bins, 
M max (logr), in the SMC and the LMC. They show that for a 
continuous power law mass function Mmax (log r) increases with 
age. This is because for a power law mass function with index 
—a, M ma x scales with the number of clusters, N, as Mmax c< 
jV 1 /'" -1 '. Assuming a constant CFR, the number of clusters 
per logarithmic unit of age, dN/d log r, scales linearly with r, 
since diV/dlogr oc rdN/dr, with d7V/dr constant. In this case 
M m ax(logr) oc r 1 ' . For a = 2 this results in a lin ear scaling 
of M max (log r) with r. lKumai. Basu & Fuiimotol d 19931) find such 
a relation in the log(Af) vs. log(r) plane of clusters in the SMC, 
LMC, M31 and M33, which t hey attribute to time-dependen t con- 
ditions for cluster formation. Elmegreen & Efremov ( 1997) were 
the first to suggest that such an increase is purely a statistical effect 
that arises from sampling cluster masses from a power law func- 
tion with index —2. This power law form was suggested for the 
mass function of cloud cores as th e result of the fractal and turbu - 
lent nature of the interstellar gas dElmegreen & F algarone 1996). 
This scale-invariant structure of the gas combined with a near con- 
stant star formation efficiency then leads to the same simple power 

law form for the mass function of clusters (Elmegreen & Efremov 

-i 

1997). For a constant star formation rate (SFR), the number of clus- 
ters is set by the time interval that is considered, which then deter- 
mines the maximum cluster mass that can (on average) be expected. 
What this means is that the probability of finding the physical con- 
ditions, such as density and pressure, to form a massive cluster is 
determined by the structure of the gas, which can be described by 
a sim ple functional form . 

IHunter et all d2003l) find M max (logr) oc r 7 in the SMC 
and LMC, from which they conc lude that a ~ 2. 4 . How ever, 
from direct fits to t he CI MF em p Ide Griis & Andersl d2006l) and 
Ide Griis & Goodwin! d2008l) find indices of a = 1.8 ± 0.1 and 
a = 2.00 ± 0.15 for the LMC and the SMC, respectively. The 
difference between the results of Hunter et al. and de Grijs and 



co-workers can be explained w hen a Schechter d istribution func- 
tion (equation [T) is assumed: Hunt er et alj d2003l) only use the 
most massive clusters, i.e. close to M* where the mass function 
is steeper, whereas the studies of de Grijs et al. use the full mass 
range of only the youngest clusters, for which Mmax < M*. From 
this it can be concluded that in these galaxies it takes more than 
10 Myr to sample enough clusters from t he CIMF to reach M t . This 
scenario is confirmed by the results of Ide Griis & Andersl d2006l) 
who show that the mass function of LMC clusters in different 
age bins gets steeper with age, from —1.8 ± 0.1 at ~ 10 Myr to 
— 2.2±0.1 at ~1 Gyr. Because these mass functions are constructed 
in equally spaced logarithmic age bins similar arguments hold as 
above: in the older age bins longer time intervals are considered, 
thus sa mpling higher masses from the CIMF where it is steeper. 
Lastly, iLarsenl d2008l) recently showed that a Schechter function 
with M* = 2 x 10 5 Mq provides a good description of the high- 
mass end of the mass function of LMC clusters when the full age 
range i s considered (but exclu ding the globular clusters). 

iGieles & Bastianl d2008l) also looked at M max (log r) for clus- 
ters in the SMC and LMC (using the data from Hunter et al.) and 
added cluster populations from five other galaxies. A linear in- 
crease of Mmax (log r) with r holds for the Milky Way open clus- 
ters and for the clusters in the SMC, LMC, M33 and M83 up to 
ages of ~ 100 Myr, which implies that a — 2 for the mass range 
(sampled in that age range) in those galaxies. It also means that 
there has been no disruption of massive clusters in these galaxies in 
this ag e range. In the 'universal' MID model of Whitmore et al. 
(2007), 90% of all clusters gets destroyed each age dex, which 



predicts a constant number of clusters in logarithmic age bins (if 
dN/dr oc r _1 , then d N/d log r=constant), an d therefore a con- 
stant Mmax in such bins dGieles & Bastianl l2008). Since this is not 
the case in most of the galaxies, the linear scaling of Mmax (log t) 
w ith r is an important arg ument against the 'universal' MID model 



w itn r is an important argu 
of IWhitmore et alj ( 120070 . 



When the full age range of the cluster population is consid- 
ered, the increase of M max (lo gr) with r is slower t han linear, 
roughly M raax (logr) oc r ' 7 dGieles & Bastianl 120081) . in agree- 
ment with what Hunter et al. had found. For clusters in M51 the 
increase of Af max (logr) with r is even slower and for the Anten- 
nae galaxies Mmax (log r) is essentially independent of r, suggest- 
ing that due to the high SFR of these galaxies Mmax — M» already 
at short intervals of r, such that A/ raax (log r) is approximately con- 
stant. However, a high MID fraction and a truncation in the CIMF 
have m ore or less the same eff ect on the evolution of Af max (log r) 
with t (Gi eles & Bastianll2008l) . In addition, a non-constant forma- 
tion history of clusters, one which was lower in the past, can also 
produce a constant Mmax (log T )- Additional age and mass distribu- 
tions are needed to tell these effects apart. 

In Section|4]it is shown that the truncated CIMF scenario com- 
bined with mass-dependent disruption nicely explains the age dis- 
tribution and mass function of M5 1 clusters and the MID model 
is rejected based on the age distribution of massive clusters. How- 
ever, a truncation in the mass function does not explain the r _1 age 
distribution of the Antennae clusters. 



2.1.2 Variations in the cluster luminosity function 

Indirect evidence for a steepening of the cluster mass function at the 
high-mass end comes from the luminosity function (LF) of clusters. 
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The LFs of clusters in different galaxies can be well approximated 
by a power law, but with an index s maller than the index of the 
CIMFemp (between -2 . 5 and -2, e.g . lDolphin & Kennicutl 20021 : 
Elme green et al.l 1 20021 ; lLar sen 2002) wi th the LF being steeper 
at higher luminosities i Whitmore et al.l 1 19991; Zepf et all 19991; 
Bened ictet alj2002l;lLarsenl2002l : lMengel et al J 20051 ; I Gieles et al.1 
2006a bc lHwang & Led2008h . The LF consists of clusters with dif- 
ferent ages. Due to the age dependent light-to-mass ratio of clus- 
ters, the LF does not necessarily have the same shape as the mass 
function. However, if the CIMF is a continuous power law with the 
same index at all ages, i.e. with the physical maximum much higher 
than A-/ max , the LF will be a power law with the exact same index. 
Age dependent extinction or bursts in the formation rate would not 
cause a difference between the CIMF and the LF. An addition of 
identical power laws always results in the same power law. So the 
fact that the LF is slightly steeper than the CIMF emp is already a 
strong indication that the CIMF is not a continuous power law func- 
tion. 

When an abrupt truncation of the CIMF at some mass M up 
is assumed, it is possible to roughly estimate the index of the 
bright-end of the LF. Assume that the CIMF is fully populated, 
i.e. the mass of the most massive cluster actually formed, Af max , is 
equal to M up . Then assume a constant formation history of clus- 
ters, so a constant number of the most massive clusters per unit 
of time: diV^/dr =constant. The luminosities of these clusters, 
L up , however, depend strongly on age. The light-to-mass ratio, 
or the flux of a cluster of constant mass, scales roughly with r 
as t _i \ with 0.7 < ( < 1 depending on the filter, such that 
dL np /dr oc 



' x oc L]^ 1 ^. Then the LF of such clusters is 



diVi 



up 



dN up d', 



dL u 



dL up 

-i/C 



(2) 



oc L~ 2 ' 5 , (4) 

where in the last step ( = 0.7 is used (BL03. lGieles et al.ll2007h . 
such that the index of —2.5 holds for the V-band LF. The same 
arguments hold for the luminosities of the 2 nd , 3 d , etc. most massive 
clusters, such that the bright-end of the LF of the entire population 
is an addition of power law with index —2.5, resulting in a power 
law with index — 2.5 . 

The faint-end of the LF should still be a power law with index 
—2 or flatter if MDD is important. This dou ble power law shape 
for the LF was found bv lGieles etafl feo06b) and I Whitmore et all 
(1999) for the LF of clusters in M51 and the Antennae, respec- 
tively. When a Schechter function for the CIMF is considered, the 
logarithmic slope of the LF shows a smooth decli ne from —2 t o 
roughly -3 between M v ~ -4 and Mv ~ -12 fcarsenll2008h . 
MDD affects the faint-end (M v > -8) of the LF and makes it 
slightly shallower (index > — 2). 



2.1.3 The value of M* 

Is the value of M* universal? Or does it somehow depend on 
galactic conditions, such as the SFR? The two proxies for the 
shape of the high-mass end of the CIMF, namely the evolution 
of M m!lx (logr) with r (Section 12.1.11 ) and the power law index 
of the LF (Section [2.1.2K both show more convincing signatures 
of a truncation in the mass function in galaxies with a high SFR. 



This suggests that in such environments there are numerous clus- 
ters with masses above M*, as opposed to galaxies with a low 
SFR. This excludes a linear dependence of M* on the SFR. If A/* 
would scale linearly with the SFR, then it would take the same 
amount of time in each galaxy for A/ max to reach A/*, since in 
the power law regime of the Schechter function M max (loRr ) oc r 



(Secti on 12. Lit. This was suggested by Weidner. Kroupa & Larsenl 
J2004I) and|Maschberger & Kroupa fc007h . who claim a formation 
epoch of lOMyr for an entire CIMF. If the CIMF is populated 
up to the highest mass each lOMyr, than there should be no in- 
crease of A/ max (log t) with age (assuming all age bins have a width 
larger than lOMyr). Bu t this is not w hat Hunter^tal] d2003l) and 
iGieles & Bastianf d2008h find. IWeidner et al.l d2004h base their re- 
sult on the relation between the luminosity of the brightest cluster 
in a galaxy, M 1 b " ghtel1 ', with the SFR of that galaxy using data from 



ILarsenl d2002h and the assumption that the brightest cluster is also 
the most massive cluster. They then use a constant light-to-mass 
ratio to convert lu minosities to m asses. This assumption was fur- 
ther investigated by ILarsenl 1 2008), who determined the ages of all 



brightest 



clusters from Larsenl d2002h and showed that the brightest 



M, 

clusters (M^ 18hle!,t < — 11) are young (<30Myr) while the clusters 
with A/^' Ehle!,t > —11 have a large age spread (between ~10Myr 
and ~ 1 Gyr). If the age of My lghit&l is My dependent, then so is 
the light-to-mass ratio. Assuming that it i s const ant can thus lead 
to a misinterpretation of the data. ILarsenl {2008) showed that the 



M- 



brightest 



vs. SFR relation for (quiescent) spiral galaxies is consis- 



tent with a constant M* of M, ~ 2 X 10° M Q . 

This value of A/* is probably not universal, since there are 
extremely massive cluste rs known (10 7 — 10 s M^. lMaraston et all 
l2004l;lBastian" et alfeooil) . which are unlikely to have forme d from 
a Schechter function with M, = 2 x 10 5 M Q . Jordan et al. ( 2007, 
hereafter J07) determine M* values of globular cluster popula- 
tions of early-type galaxies in the Virgo Cluster by comparing 
the cluster luminosity functions to 'evolved Schechter functions'. 
They find a trend of M» with host galaxy luminosity, with Af» 
being higher in brighter galaxies, and values for M, between a 
few times 10 5 Mq and a few times 10 6 Mq. Several other studies 
have also shown that a Schechter function fits the high-mass end of 



law with index —2 i 


McLaughlin & Pudritzlll996l;lBurkert & Smith! 


l200d;lFall& Zhang! 


2001 ; Waters et alj2006l; Harris et alj2008l). 



iBillett et"al] d2002h s uggest the r e shou ld exist a maximum 
cluster mass based on the Kennicutt| d 19981) relation between the 
SFR and gas density and the assumption of a pressure equilibrium 
between the ambient interstellar medium and the cloud-cores from 
which star clusters form. For a constant cluster density this rela- 
tion is Mup oc SFR 2 and for a constant cluster radius, it would be 
M up oc SFR 2/3 . With the latter increase of M up with the SFR, i.e. 
slower than linear, it is possible to explain the observational trends 
discussed in Sections [2. 1 .1442. \2\ The key point is that M* is prob- 
ably environment dependent and for the masses of a cluster popula- 
tion in a given galaxy to be affected by the truncation, the product of 
the SFR times the age range has to be high. For A/* oc SFR 2 / 3 , low 
SFR galaxies need longer time intervals for M max to reach M* than 
high SFR galaxies, which seems to be what is found from observa- 
tions. The relation between Af max and the SFR will be quantified in 
more detail in Section|2~2l 
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T-SFR-At [M ] 

Figure 1. Example of the relation between the mass of the most mas- 
sive cluster, Af max , as a function of the total mass formed in clusters, 
T ■ SFR ■ At, when the CIMF is a Schechter function (equation [TJ with 
M, = 2x 10 s M Q lLarsenll2008l) . The dashed line shows a prediction 
for A/max in the case of an nontruncated mass function. The dotted lines 
indicate three masses and the corresponding logarithmic slope, —j3, of the 
Schechter function. The top x-axis shows the time range that is needed to 
form the amount of mass that is represented on the bottom rr-axis when 
SFR = 1 AZoyi-- 1 and V = 0.1. 



2.2 The relation between A/ max and the star formation rate 

The mass of the most massive cluster in a cluster population, M max , 
depends on how many cluster are formed. This can be related to 
the CFR which scales linearly with the constant A in equation QJ, 
since 



CFR = 



f 00 dN 
Jm_ dMi 



AEi 
10 A, 



Mwm 



(5) 

(6) 
(7) 



shape of the age distribution, which can be derived once the CIMF 
is evolved with disruption (Section^. 

The relation between A/max and A can be found from 



1 = 



dN_ 

dMl 



AE 2 (M ma x/M, 
Afmax 



(8) 



(9) 



For M m ax << M* the term E 2 (M max /M*) ~ 1, and A ~ M max , 
which is the result for a continuous power law with index — 2. An 
expression for A/ max as a function of the CFR can be found by sub- 
stituting A as a function of the CFR from equation l[6ll in equa- 
tion (|9). Since this is then in units of AfQyr -1 , a multiplication by 
some time interval, At, is necessary to get an expression for Af max 
in Mq. Assuming that a constant fraction F o f the SFR end s up 
in star clusters that survive the embedded phase (Bastian 2008j), i.e. 
CFR = r ■ SFR, then a relation between M max and the SFR is found 



A/ m 



T ■ SFR ■ At 

Ei (M mm /M«) 



(10) 



E 2 (M max /M < ,) 

This relation between M ma x and T ■ SFR • At is illustrated in 
Fig-H Since the exponential integral is not easily inverted the prod- 
uct T • SFR • At is calculated for a series of M max values. A value 
of M, = 2 x 10 5 M (7) is used, which was found for spiral galax- 
ies bv lLarserj (2008), and Ei(A/ m i„/M*) = 10, corresponding to 
Mmin ~ 1. The dotted lines indicate three masses where the loga- 
rithmic slope of the mass function, — j3, has the values —2.1,-2.3 
and —3. The logarithmic slope at M of a mass function dTV/dAf is 
defined as 



dln(diV/dM) 



dlnM 
which for the CIMF gives 

— P(Mi) = — 2 — Mi/M*. 



(11) 



(12) 



Note that j3 is defined such that it is positive for a declining 
dN/dM. 

From Fig.Q]we see that M max ~ 0.1 V ■ SFR • At for A/ max < 
0.1 M*. This corresponds to the mass range of the CIMF where 
f3 ~ 2. For higher values of V ■ SFR • At the relation between 
A/max and T ■ SFR • At flattens. The dashed line shows the predicted 
Afmax for a nontruncated mass function, or a much higher M» . Even 
though the two predictions diverge to a difference up to two dex in 
A/max, a large increase of T ■ SFR ■ At (three dex) is needed to reach 



is found, showing that the relation between the CFR and A is rel- 
atively insensitive to the ratio Af ra i n /M». The variable A can be a 
function of time, A(t), to include variations in the CFR. This does 
not influence the shape of the mass function, but it will affect the 



where E n (x) with n = 1 is a generalised expression of the expo- this difference. When M max in the CIMF cmp is sampled close to or a 
nential integral. In the last step M mm /M, ~ 10" 5 is used. For a bit above M * > the difference between a continuous power law and a 
ratio of M mi „/Af, = lO^lO" 3 ) a relation of CFR ~ 8A(6A) Schechter function will be hard to tell. Mainly because the number 

of clusters in this high-mass tail is low. 

The top x-axis of Fig. Q] is labelled with At values corre- 
sponding to r ■ SFR ■ At values with a fixed SFR and V of 
SFR = lMQyr - and T = 0.1. These At values can be in- 
terpreted as age ranges. In this illustrative example of a galaxy 
with a moderate SFR, the CIMF erap (At < lOMyr) will con- 
tain clusters with masses up to ~5xlO 4 M0, while clusters with 
masses up to ~10 6 Mq have formed when an age range of 
~10Gyr is considered. However, without the exponential trun- 
cation clusters with masses up to 10 s Mq are predicted to have 
formed. The example in Fig. Q] can be applied to the Milky 
Way. Its most massive young star cluster known to date is West- 



2 The exponential integral is defined as E n = t~ n exp(— xt)dt, and 
is related to the (upper) incomplete Gamma function, T(a, x), as E n (x) = 
x 71-1 r(l — n, x). See Section 6.3 in Numerical Recipes iPress et al] 
[1993) for details on the exponential integral and how to implement this 
function in a code. The exponential integral is predefined in IDL as the 
function expint (n, x) . 
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Figure 2. Different stages of the evolved Schechter mass function for a constant mass loss rate (left panel, equation 1171 ) and a mass-dependent mass loss rate 
(right panel, equation 1261 ). The corresponding values of the ratio A/M* (left) and A 7 /Af^ (right) are indicated. The dashed line in both panels shows the 
Schechter function at t = (equation[T}. 



erlund 1 which has a mass arou nd 10 5 M© dClark et all 120051 : 
iMengel & Taccon i-Garman 120071) . There are a handful of other 
massive clusters in the Milky Way known, such as the Arches and 
Quintuplet cluster towards the Galactic centre, NGC 3603 and the 
two recently disco vered red super giant clusters (Fig er et ai]|2006l ; 
iDavies et alj2007h . They are all around ~10 4 Mq and span an age 
range of ~10 Myr. From Fig.Q]it can be seen that for this age range 
values of 1 — 5 x 10 4 M@ are expected for A/ max . If the last Gyr 
of cluster formation is considered, the Galaxy has probably formed 
several clusters more massive than Westerlund 1, but not as many, 
nor as massive compared to the scenario in which the CIMF is a 
continuous power law. 

With the properties of the CIMF introduced, an analytical ex- 
pression for the 'evolved Schechter function', i.e. the cluster mass 
function after mass-dependent disruption is applied, can now be 
derived. 



3 THE EVOLVED SCHECHTER FUNCTION 

In this section the 'evolved Schechter function' is presented, based 
on the CIMF of equation {T} which is evolved with mass-dependent 
cluster disruption due to e.g. the tidal field of the host galaxy and/or 
encounters with GMCs. Cluster mass loss due to stellar evolution 
is not considered at this stage. Also, the power law part of the mass 
function is always assumed to have an index of — 2 for the moment. 

It is a simple mathematical exercise to derive all expressions 
in this section for a variable index and to include the effect of mass 
loss by stellar evolution, but for reasons of simplicity this is not 
done in the derivations of the formulae in the coming two sections. 
In Section[3~3l these two effects are added to the formalism. 



constant. The disruption time for a cluster of mass M, i.e. the time 
needed to reach M = 0, is then given by 



tdis 



M 

to M, 



(13) 
(14) 



where in the second step to = 1/ M is introduced. This linear scal- 
ing of tdis with M follows from the assumption that tdis is a con- 
stant times the half-mass relaxation time, t T h, and the assumption of 
a constant cluste r half -mass density and a constant Coulomb lo ga- 
rithm in t rh (e.g. iGnedin & Ostriker|[l997l ; iFall & Zhang|l200lh . A 
more general expression for the scaling of tdis with M of equa- 
tion d 14b would be tdis = toM 1 (equationl22t. with 7=1. 
The mass evolution of a cluster with tim^fl is then 



M 



M, - Mt, 
Mi - A, 



(15) 



where the variable A = t/to(= Mt) is introduced, which is a 
proxy of time, but has the dimension of mass since it is the amount 
of mass with which Mi has reduced at time t due to disruption. 

Though the assumptions that have to be made to arrive to this 
constant mass loss are arguable, the result is a mathematically ap- 
pealing description of the evolution of cluster masses, enabling us 
to make simple analytical predictions. In Section [3~2l a more realis- 
tic M, based on results of iV-body simulations, is explored. 

The mass function as a function of time, diV/dM, follows 
from the conservation of number 



diV 
dM 



d7V 
dMi 



8Mr 



dM 



(16) 



3.1 Constant mass loss rate 

The evolved Schechter function is first introduced in the form de- 
scribed by J07. The mass loss rate of clusters, M, is assumed to be 



3 To keep the formulae synoptic, the variables M and AN/ AM, i.e. without 
the subscript i, are used for to the mass and the mass function as a function 
of time, i.e. M(t) and AN/AM(t), respectively. The variables M and M% 
are in units of Mq and t and to are in units of time (Myr). 
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Since dMi/dM 
follows that 



1 and Mi = M + A (equation [BJl it then 



diV 
AM 



A 



[M + A] ' 



• exp 



M + A 

~mI 



(17) 



The behaviour of this evolved Schechter function is shown in the 
left panel of Fig. [2] The dashed line shows the CIMF (equation (T) 
or equation 1171 with A = 0) and the full lines show the result of 
equation (TTJ at different times (i.e. for various ratios of A/M*). 
For M < 0.1 A and A < M* the AN /A M is flat (/3 = 0). T his is 
the direct consequence of the constant M dFall & Zhand200ll) . The 
reader is referred to J07 for an equivalent functional form describ- 
ing the luminosity functions, i.e. the number of clusters in constant 
magnitude interval, which is peaked and rises at low luminosities. 

From equation d!7t an expression for the logarithmic slope 
(equation II It can be derived 



■f3 = -2 + 



2A M 
M + A ~ MZ 



(18) 



The three terms on the right-hand side are, respectively, the initial 
index of the power law part of the Schechter function, a disruption 
term that makes the mass function shallower and a truncation term 
that makes the mass function steeper. Solving for M gives the mass 
at which the mass function has the logarithmic slope —/3 



M(/3) = . 



-A + (P-2)M t + ^/(A+{f3 - 2)M*) 2 +8AM, 



(19) 



From this it can be seen that M (ft) assumes its simplest form when 
P = 2. This is where the truncation term in equation J 1 3b equals the 
disruption term and this could, therefore, be interpreted as a point 
in the mass function where there is a balance between a steepening 
due to the truncation and a flattening by disruption. Using equa- 
tion l |19t with f3 = 2 an expression for M(/3 = 2) arises 



M(/3 = 2) = 



-A + VA 2 + 8 AM. 



(20) 



From equation d20b it can be seen that in the limit of A << M», 
i.e. at young ages where (massive) clusters are not yet affected by 
disruption, the scaling of M(/3 = 2) with A converges to M(/3 = 
2) = %/2AM*. Since A is a proxy for age, it means that the part 
of the mass function that can be approximated by a power law with 
index —2 is found at higher masses at older ages. This can also 
be seen from the expression for (3 as a function of M and A in 
equation dl8t . For the mass function to have a logarithmic slope of 
P — 2, the two right-hand terms (i.e. the disruption and truncation 
term) in equation d 1 81 ) need to have a sum of 0, so A is higher at 
higher masses. 

From equation l |19t it can be seen that in the same regime of 
A << M» the behaviour of M(f3 / 2) is different. For j3 < 2 , 
M(/3) also increases with A, but linearly, so faster than Mifi = 2) . 
For (3 > 2, M{f3) is constant. This is because at young ages the 
steep part of the mass function is not affected by disruption at all. 

Equation dl9t also allows us to find the mass where the mass 
function, when presented as the number of clusters per logarithmic 
unit of mass, AN/ A log M, peaks or 'turns over' . This is where /3 = 
1 and this M is defined as the turnover mass, Mro- Its dependence 
on aforementioned variables is (J07) 



Mr 



-(M, + A) + ^/(A + M t ) 2 + 4AM* 



For A << M* the turnover mass scales as Mto oc A, like M(/3 < 
2), meaning that the turnover shifts to higher masses as time pro- 
gresses. Notice that due to the different scaling of M(/3 = 2) and 
Mto with A, they approach each other for A/M* > 1. In fact, 
for A >> M* the shape of the mass function remains unchanged 
(J07) and approaches an equilibrium form with Mto = M» and 
M(/3 = 2) = 2 M». The number of clusters reduces as time pro- 
gresses, but the shape of the evolved Schechter function remains the 
same. This means that clusters with initial masses much higher than 
M* are replacing clusters with lower masses that have already been 
disrupted. Since the number of clusters above M, is very low, in 
practise it means that the total number of remaining clusters quickly 
drops to after A = M* . 

The behaviour of the relations described above is illustrated in 
the left panel of Fig. [3] It shows M vs. A, with both quantities rel- 
ative to M* , such that the x and y scales are in dimensionless units 
of mass and time, respectively. The results for/3 = 1/1.7/2/2.3 
(equations (19), f20) & Ell ) are shown. 

The grey shaded region in Fig. [3] shows that there is a large 
region in the age vs. mass plane, where the logarithmic slope of 
the mass function equals — 2 ± 0.3. This region more or less coin- 
cides with the strip that observed cluster populations occupy in the 
age-mass plane, implying that the observable part of the evolved 
Schechter function can be approximated by a power law with index 
—2. This will be verified in Section [4] using empirically derived 
ages and masses of star clusters in M51. But first the case of a 
mass-dependent mass loss rate is discussed in Section l3~2l 



3.2 Mass-dependent mass loss rate 

In the previous section a (simplified) linear scaling between idis 
and tib was adopted. However, there are theoretical arguments for 
a de- coupling of idis from i r h for clusters dissolv ing in a tidal 
field dFukushige & Heggiel2 000). Bau mgardtl d200ll) predicted that 
tdis oc 4 for clusters that are initially Roche-lob e filling. This 
leads t o a simple scaling of idis with M of the form dLamers et al.l 
2005b) 



idis =t M 7 , 



(22) 



with to an environment dependent constant and j ~ 0.65. 
This result for Roche-lobe filling clusters was indeed found 
from iV-body simulation s of c l usters dissolving in a tida l 
field dVesperini & Heggid 1 19971 ; iBaumgardt & Makinol 120031) 
and it was also found for Roche-lobe underfilling clusters 
dGieles & Baumgardtl 120081) . The simple disruption law of equa- 
tion \22i was also used to describe age and mass distributions of 
clusters, from which a mean value of 7 = 0.62 ± 0.06 was found 
(BL03). 

The mass loss rate, M, can be related to idis as 



M 

M = — , 

Idis 

M 



l- 7 



to 



(23) 



(21) 



For 7 / 1 this results in a mass loss rate that is no longer con- 
stant, but instead becomes a function of M. Note that then to is 
not simply 1 / M anymore and that idis is no longer the time of total 
disruption as in Section [3~T| The total disruption time, i.e. the time 
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Figure 3. The evolution with age of the mass where the evolved Schechter function has the logarithmic slope, —j3, for different values of (3. In the left panel 
a constant mass loss rate is considered (equation [T7] Section |3Tt and in the right panel the mass loss rate depends on mass (equation |26| Section [T5J. The 
a; -axes represent time or age since A = t/trj and A 7 = ■yt/to. The mass for which (3 = 1, i.e. the mass at which the mass function, presented as the number 
of clusters in logarithmic mass bins, turns over, is indicated as M-ro- The grey shaded region indicates the part of the mass function where the logarithmic 
slope is -2.0 ± 0.3. 



it takes for all stars to become unbound, will be referred to as t^' s . 
The relation between t'^ s and t<jis will be derived a bit further down. 

With tdis denned as in equation d22t it is still possible to 
jet an analytical ex pression for the mass evolution with time 



get an analytical exp 
dLamers et al1r20 05a) 



M = [M 7 - A. 



iV-y 



(24) 



where A 7 is defined as A 7 = ft /to to get a similar looking ex- 
pression as in the constant mass loss case (equation [T3J- The effect 
of stellar evolution can be taken into account by replacing M 7 in 
equation d24t by (/JbM,) 1 , where /i cv is a time-dependent vari- 
able that represents the fraction of the initial stellar mass that has 
not been lost by stellar winds or super-nova explosions (typically 
between 0.7 and 1). For reasons of simplicity ^t C v is not consid- 
ered in the derivations that follow in this section. In Section [3~3l the 
results including /i ev are given. The symbol A 7 is used when refer- 
ring to 7 < 1 and A is used for the case of 7 = 1. The variable A 7 
is a proxy of time, like A, but note that A 7 does not represent the 
amount of mass that is lost from Mi, like A, because of the powers 
of 7 and I/7 in equation d24t . The instantaneous disruption time 
for a cluster with mass M is when A 7 /M 7 = 1 (equation 124b. 
such that = (to/7) M\ or t m = 7 C- 

The partial derivative that is needed (equation 1 161) to get an 
expression for the evolved Schechter function can be found from 
equation d24t 



dMj 
~dM 



1/7-1 



(25) 



which combined with equation (TJ & d!6t gives an expression for 
the evolved Schechter function for a variable 7 



diV 



AM" 1-1 



AM [Mt + A. 



l(7+l)/7 



exp 



[M 7 + A 7 ] 1/7 



Ah 



(26) 



For M « A'h the exponential term in equation i26\ is about 1 
and the result for an nontruncated mass function of lLamers et al.l 
( I2005an is found, which at low masses and old ages (M 1 « A 7 ) 
goes to AN I AM oc M T_ . This is the consequence of the assumed 
relation between M and tdis (equationl23t. If M was assumed to be 
constant in time, even when ^(Mj) oc M 7 , with 7 < 1, the low- 
mass end of the mass function would be flat (index of 0). The shape 
of the low-mass end of the evolved mass function is thus indepen- 
dent of the shape of the CIMF a t those masses and depends only 
on the way clusters lose mass jFall&Zhangll200ll ;l Lamers et al.l 
l2005ah . 

Using the definition of the logarithmic slope, — /3, from equa- 
tion (TTJ an expression for /3 is found 

( 7 +l)A 7 M 7 (M 7 + A 7 ) 1//7_1 



(27) 



Mi + A 7 M, 

Unfortunately equation d27b can only be solved for M analyti- 
cally for a few specific values of 7 and /3 and the solutions are rather 
complicated. Therefore, M(f3) is solved numerically and the right 
panel of Fig. [3] shows the result for the same values of f3 as in the 
constant mass loss case (left panel of Fig. [3j( using 7 = 0.65. For 
M << M„, Mto scales as Ay 7 , which was also found under the 
assumption of a continuous power law in the instantaneous disrup- 
tion model of BL03. However, the increase of Mto with age slows 
down at high ages, due to the truncation. In contrast to the constant 
mass loss case, the evolved Schechter function with 7 < 1 does not 
approach an equilibrium shape. For 7 = 0.65, the value of Mto 



reaches its highest value (Mto — 0.2 M») around A 7 



2M? 



and then decreases again for older ages. This subtle difference 
might be important. If globular clusters formed from a Schechter 
type CIMF, with M* ~ 10 6 M , then together with the fact that 
Mto — 2 x 10 M©, the constant mass loss model would say that 
the value of Mto can still increase by roughly a factor of five (see 
left panel of Fig. [3}- However, the 7 = 0.65 model shows that Mto 
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already has its highest possible value when Mjo/M* ~ 0.2 (right 
panel of Fig.O. A possible explanation for the near universality of 
the value of Afro could thus be that most globular cluster systems 
have reached already their maximum value of Mto /M, due to dis- 
ruption and for a range of roughly an order of magnitude in A 7 , the 
value of Mto will be (within a factor of two) about 0.2 M*. Since 
globular clusters are roughly coeval, a spread in A 7 can be inter- 
preted as a spread in disruption time-scales, thereby allowing for a 
range of tdis values (i.e. different galactocentric distances) resulting 
in similar values of Afro - Perhaps this adds a piece to the puzzle of 
the probl em of the universality of the globular cluster mass function 
(see alsolOstriker & G nedin 1997; Vesperini 200d: IWhitmore et al.l 
l2002l ; lBaumgardt et al.ll2008l) . 

The mass for which the mass function has an index of —2 
scales approximately as A2, with rj ~ 0.6, i.e. slightly different 
than the vA scaling found for the 7 = 1 case described in Sec- 
tion [3J0. However, the A 7 6 scaling is a bit closer to the increase 
with age of the limiting mass due to the detection limit, Miim(r), of 
empirically derived cluster masses. For a V-band detection limit, 
Miin,(r) scales as r ' 7 . This implies that the observable part of the 
evolved Schechter function, i.e. above the detection limit, has a log- 
arithmic slope of approximately —2 at all ages. This is an important 
result, because it means that there is no significant difference to be 
expected between the CIMF emp and the shape of the cluster mass 
function at old ages. So the argument that disruption needs to be 
mass-independent because old clusters have a similar mass func- 
tion as young clusters does not hold when a Schechter function is 
assumed for the CIMF. 

All results presented for the constant mass loss case in Sec- 
tion |3.1| are simply a subset of the more general solutions presented 
in this section. That is, when using 7 = 1, A 7 = A and equa- 
tions J22l >. l !24t . d26t & ( 1271 ) from this section are the same as equa- 
tions fl4l >. < fl~5T >. JTTt & (Q~8}, respectively, from Section 13.11 For 
completeness a general expression for the evolved Schechter func- 
tion, including mass loss by stellar evolution and a variable initial 
power law index, is presented in Section [3~3"1 



( 7 + a-l)A 7 M 7 (M 7 + A 7 ) 



1/7-1 



(29) 



M~< + A 7 flevM* 

From equations d28t & ( 1291 ) it can be seen that /i ev does not 
affect the shape of the mass function at low masses, it only affects 
that vertical offset by a bit. The exponential truncation occurs at 
slightly lower masses, at (i„M, instead of M*. 



4 APPLICATION TO M51 

To illustrate the cluster population model of the previous section, 
some of its elements are compared to the cluster population of 
the interacting galaxy M51 (NGC 5194). The cluster ages and 
masses were determined from HST/WFPC2 multi-band photom- 
etry by Bastian et al. (2005, hereafter B05). The masses are cor- 
rected for mass loss by stellar evolution, i.e. they represent the sum 
of the initial masses of the stars that are still bound in the clus- 
ters. The present day masses of the oldest clusters are, therefore, 
~25% lower. The M5 1 cluster population is a good benchmark to 
test the evolved Schechter function for several reasons: firstl y, it is a 
rich c luster population host i ng several massive st ar clusters dLarsenl 
l2000l : B05: lLee et ai1l2005l ; lHwang & Leell2008l) . M51 hosts suffi- 
cient clusters above ~10 4 Mq, which is the approximate minimum 
mass for the optical fluxes not to be affected by stochastic effects 
due to sampling of the stellar IMF, allowing reliable age datin g 
through multi-band photometry (e.g. ICervi no & Luri dianj|2004t) . 
Secondly, for this cluster p opulation it has been suggested that th e 
mass function is truncated iGieles et al.ll2006bl : lHaas et al.| [2008). 
Lastly, the di sruption time of clu sters in M5 1 is thought to be very 
short (BL03. iGieles et al.ll2005l) . due to the strong tidal field and 
the high density of GMCs in the inner region of the galaxy, where 
most of the clusters of the B05 sample reside. These environmen- 
tal properties suggest a noticeable evolution of the mass function. 
Before a comparison between the model and the data is given, the 
parameters for the CIMF (M») and the disruption law (equationl22t 
are determined in Section|4. II 



3.3 Including stellar evolution and a variable power law 
index 

As mentioned in the beginning of this section all derivations are 
done for a fixed power law index of —2 (equationQJ and without the 
inclusion of mass loss due to stellar evolution. The expressions for 
the evolved Schechter function and the logarithmic slope —(3 are 
here given for a variable initial power law index, —a, and an addi- 
tional term which is the fraction of the original mass that is not 
lost due to stellar evolution. Therefore, is a function of time and 
can be taken from an SSP model. See also lLamers et al.l fc005al) for 
analytical approximations of /i cv (t) . Following the same steps as in 
Section [372"l the expressions for dTV/dAf and (3 are 



dTV 
dM 

and 



7-1 



[Mt + A. 



l( 7 + Qf-l)/ 7 



exp 



[AP+ A 7 ] 



1/7 



(28) 



4 In fact, by numerically solving equation {27} for /3 = 2 and < 7 1 
the relation between 77 and 7 is r\ ~ 1 — 0.85 7 + 0.35 7 2 . 



4.1 A fit of the evolved Schechter function to the M51 cluster 
population 

If a constant cluster formation history is assumed the evolved 
Schechter function (equation I26t essentially represents the prob- 
ability of finding a cluster with mass between M and M + dM and 
age between r and r + dr. It can thus be used as a two-dimensional 
distribution function to determine all the parameters of the evolved 
Schechter function from empirically derived ages and masses us- 
ing a maximum likelihood estimate. The only thing which needs 
to be added is the detection limit. For this the result of the in- 
completeness analysis of B05 is used, who show that their cluster 
sample is limited by a detection in the F435W (roughly B) band, 
for which the 90% completeness fraction is a 22.6 mag. The same 
simple stellar population (SSP) model was used as in B05 to de- 
rive the ages and mass, namely the GALEV models for Salpeter 
stellar IMF between 0.15M B and 50 M B dSchulz et alj I2OO2I; 
lAnders & Fritze-v^ Alvensleben||2003l) . The photometric evolution 
of a single mass cluster in the F435W band from the SSP model, 
combined wit h the 90% completene ss limit and the distance to 
M51 (8.4Mpc. lFeldmeier et alj|l997h then giv es the limiting clus - 
ter mass as a function of age, A/ij m (r) (BL03; lGieles et al]|2007l) . 
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This limiting mass, together with the cluster ages and masses, is 
shown in Fig.[4]and in panel (a) of Fig. [5] 

Using Afii m (r) and equation l !26t artificial models for varying 
7, to and Af* are built and a (simultaneous) maximisation of the 
likelihood of these parameters gives 

• M» = (1.86 ±0.52) x 1O 5 M ; 

• t = 0.19±0.10Myr; 

• 7 = 0.67 ±0.06. 

The (statistical) uncertainties on each variable are determined us- 
ing a bootstrap method. For this the original data-set is randomly 
re-sampled 1000 times, allowing for multiple entries of the same 
age-mass pair and omission of values, such that the total number of 
age-mass pairs is the same. The likelihood of the three parameters 
is determined for each of these 1000 samples and the standard de- 
viation of the results for each parameter is used as the uncertainty. 

The disruption parameters (to and 7) imply a total dissolution 
time for a cluster with an initial mass M* of t'°' = (to / 7) M2 — 
950 Myr, implying that there should not be many clusters that sur- 
vive longer than a Gyr. Indeed B05 found only a handful of clus- 
ters older than 1 Gyr, although this is probably not only caused 
by disruption, but also by detection incompleteness. A cluster with 
an initial mass of 10 4 M© gets completely destroyed in tj' ~ 
130 My r. This latter value is in perfect agreement with the re- 
sults of lGielesetall d2005f) . who found 100 < tf < 200, de- 
pending on what is assumed for the cluster formation history. This 
very short lifetime of clusters in M51, as compared to the lifetime 
of clusters in the solar neighbourhood or the Magellanic Clouds, 
was attributed to th e high molecular cloud density in this galaxy 
dGieles et alj2006ch . 

The value of 7 = 0.67 ± 0.06 agrees perfectly with the mean 
result of BL03 who found 7 = 0.62 ± 0.06 from fits to the age 
and mass distributions of clusters in several galaxies. This value 
also follows from theory and iV-body simulations of clusters dis- 
solving in a tidal field, since for these clusters th e disruption time- 
scale can be expressed as tdis oc (N/ log A) 0,75 t Baumgard |200ll ; 
iBaumgardt & Makind 120031 : iGieles & Baumgardtl 120081) . In this 
equation, log A ~ log(O.lA) is the Coulomb logarithm and it fol- 
lows that this expression for tdis can be well approximat ed by tdis oc 
N ' 62 when a relevant range of N is considered dLamers et al.l 
2005b). The disruption time due to external perturbations scales (on 
average) in a similar way with mass. It actually scales linearly with 
the density of clusters, but since the radius of clusters depends only 
weakly on their mass, tdis due to external perturbations (GMCs, 
sp iral arms, etc.) also scales as M 7 with 7 s lightly smaller than 
1 dGieles et alJl2006d : lLamers & Gieles|[2006l) . IGieles et al.l d2005h 
also looked at the mass dependence of cluster disruption in M51 
and found that the value of 7 depends on the mass range consid- 
ered. When excluding the most massive clusters (A/ > 10 s M ) a 
value of 7 = 0.6 was found. Here it is shown that, when consider- 
ing a Schechter function for the CIMF, the full cluster population 
can be fit with a disruption law in which tdis depends on mass. 

The valu e of M* is in excellent agreement with the result of 
lLarsed d2008l) who fo und M, = (2 ± 0. 5) x 10 5 M Q for a sam- 
ple of spiral galaxies. IGieles et a l. (2006b) found a slightly lower 
value of M* = 10 s M through a model comparison to the LF of 
M5 1 clusters. Since the LF is usually constructed for apparent lu- 
minosities, i.e. not corrected for local extinction, the derived value 
of M* from the LF can be underestimating the true value of A/*. 



IGieles et al.l d2006bl) applied a correction of Av = 0.25 mag to all 
clusters to roughly account for extinction, but B05 showed that the 
youngest clusters, which are the brightest ones, are slightly more 
extincted than this. Extinction of the brightest clusters has a similar 
effect on the LF as lowering M* , s o this could be why the value for 
M* found bv lGieles et al.l ([2006b ) is slightly lower than that found 
here. 

4.2 Logarithmic slopes at different ages 

Fig. [4] shows the M51 ages and masses in dimensionless units on 
top of the predictions for the evolution of several values of j3. The 
results of the maximum likelihood estimation from Section |4~Tl are 
used for the calculation of M(/3) and the normalisation of age and 
mass. 

The increase of the minimum observable cluster mass, 
Miim(r), is shown as a dot-dashed line marking the lower enve- 
lope of the data points. Most of the clusters fall in the grey region, 
for which the predicted logarithmic slope of the mass function is 
—2 ± 0.3, with the M(/3 = 2) line showing a very similar increase 
with age as the data points. At old ages M\\ m (r) approaches the 
M(/3 = 2) line, implying that the mass function at these ages is 
somewhat steeper than a —2 power law. Although the disruption 
time of clusters in M5 1 is very short, A/to does not get closer than 
roughly one mass dex below Mu m (r) and is, therefore, not observ- 
able at any age. The dashed line in Fig.|4]shows Af (/3 = 1.7) for 
the case of a continuous power law CIMF (BL03). This line shows 
that disruption would clearly make the mass function of the old- 
est half of the data flatter than a —2 power law if the CIMF was a 
continuous power law with that index, whereas the M{fi = 1.7) 
line for the evolved Schechter function bends down at old ages and 
stays (almost) below the detection limit. 

A direct comparison between the evolved Schechter function 
(Section [3.2t and the M51 cluster mass functions at different ages 
is given in Section l4~3l 

4.3 Mass functions at different ages 

The presentation of cluster ages and masses in dimensionless 
units in the previous section is not very common. Therefore, the 
age-mass diagram is presented again, but now in physical units 
(panel [a] of Fig. |5). Three age bins are created, with bound- 
aries 4, 10, 100 Myr and 600 Myr. The upper boundaries are indi- 
cated by T\, T2 and T3, respectively. Due to the increasing A/ii m (r) 
with r, three corresponding mass limits are defined as Mj = 
Afiim(rj), with 1 < j < 3. The values of Mj are roughly 
2 x 10 3 M Q ,10 4 M Q and 6 x 1O 4 M . 

The empirical mass functions in different age bins are shown 
in panels (b)-(d) of Fig.[5] In each panel an arrow denotes the lower 
mass limit of each sample, Mj. The first mass bin starts at this 
lower limit and the number of clusters in each mass bin is counted 
and then divided by the width of the bin, such that AN /AM is ob- 
tained. The AN /AM is also divided by the age range of the sample. 
In this way the histogram points represent the number of clusters 
per unit of mass and time/age, or AN /{AM At), which can be com- 
pared to the evolved Schechter function. The variables for 7, to and 
Af* as determined in Section l4~T1 are used for the model predictions. 

The dashed lines show the CIMFs (equation[TJ and the evolved 
Schechter functions (full lines, equation [26) are based on the mean 
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Figure 4. The time evolution of the mass where the evolved Schechter func- 
tion has the logarithmic slope, —/3, for fi = 1, 1.7, 2 and 2.3 (same as in the 
right panel of Fig. [3}- The dots are ages and masses of M51 clusters from 
B05, with the dot-dashed line defining the limiting mass due to a detection 
limit of 22.6 mag in the F435W filter. The dashed line shows M(/3 = 1.7) 
when a continuous power law CIMF is assumed. The values of 7, to and 
M* , used to normalise the ages and masses to dimensionless units, follow 
from a maximum likelihood estimation (Section [4. i\ . 



age of the M5 1 clusters in each age bin, being approximately 5 Myr, 
50 Myr and 250 Myr for the panels (b), (c) and (d), respectively. 
The vertical offset is determined by A, and the best agreement is 
found for A — 0.07 M0yr -1 . The coincidence of the empirical 
mass functions with the CIMF (dashed line) in the first two age 
bins, shows that the number of clusters per linear unit of time is ap- 
proximately constant in the first 100 Myr. The implications of this 
will be discussed in more detail in Section [44] Only the dN/dM 
in the oldest age bin falls visibly below the CIMF. If a continuous 
power law distribution function would have been assumed for the 
CIMF, many more high-mass clusters would have to have been de- 
stroyed to make the model prediction agree with the data (already 
noted by Gieles et al. [ 2005] who do not include the most massive 
clusters when determining the disruption time of M51 clusters). 
The evolved Schechter functions describe the empirically derived 
mass functions very well and as could be seen already in Fig.[4] the 
turnover of the evolved Schechter functions remains well below the 
detection limit. 

The empirical mass functions are also approximated by power 
laws and their indices, —/3m, are indicated in each panel. The em- 
pirical mass functions in the first two age bins can be approximated 
by a power law with index of roughly —2, but the oldest mass func- 
tion is steeper. From a comparison to the (evolved) Schechter func- 
tion, we see this is because the limiting mass approaches M» at this 
age, where (3 — 3. On the low-mass end of the evolved Schechter 
function the effect of mass-dependent disruption is seen at old ages 
(panel [d]), where the mass function flattens to a logarithmic slope 
of 7 - 1 = -0.33 and the turnover is at M TO ^ 10 4 M . 

Though this is only an application to one data-set, it shows 
that a Schechter function evolved with mass-dependent disruption 
nicely describes empirically derived age and mass distributions. 



The lack of a flattening of the empirical mass function at old ages, 
as compared to the CIMF cmp , is due the exponential truncation at 
high masses and the detection limit that shifts the mean observed 
masses closer to M„, where the CIMF is steeper. 

4.4 Age distributions for different mass cuts 

Assuming a constant formation history of clusters, the cluster age 
distribution, dN/dr, can be found from a numerical integration of 
the evolved Schechter function from Mj to 00 at each r. 

In panel (e) three (mass limited) age distributions are shown, 
for the three maximum ages, Tj, and corresponding lower mass 
limit, Mj . The clusters are binned, such that the first bin starts at 
4 Myr, the age of the youngest cluster, and the last bin ends at Tj . 
The number of bins is chosen such that the bin widths are approx- 
imately constant for the three mass cuts. The number of clusters in 
each bin is counted and divided by the bin width, such that dN/dr 
follows, i.e. the number of clusters per Myr. 

Over-plotted as full lines are the age distribution resulting 
from the numerical integrations of the evolved Schechter functions 
again using the best-fitting parameters from the maximum likeli- 
hood estimate of Section |4~T] A single value of A — 0.07 A/oyr -1 , 
i.e. the same value as in Section l4~3l was used to describe all three 
age distributions. The different vertical offsets are caused by the 
different values of Mj. The good agreement between the model 
and the three empirical age distribution shows that both the disrup- 
tion parameters and the CIMF shape, with the parameters from the 
model fit of Section |4~T1 provide a good description of the cluster 
population of M5 1 . With the value of A and the global SFR of M5 1 
it is possible to derive T, i.e. the fraction of star formation occurring 
in star cluster s that survive the embedded phase (CFR = F ■ SFR, 
| Bastiadl2008h. The global SFR of M51 is around 5 Afoyr" 1 (e.g. 
Iscoville etai1l200lh . Since T ■ SFR = AE 1 (M min /M,) (equa- 
tion^, with Ei(Af ra i„/Af,) roughly between 6 and 10, depend- 
ing on Af ra i„, the value of T resulting from this comparison is then 
T ~ 0.08 - 0.14. So roughly 10% of the star formation in M51 
occurs in star clusters that are detectable in the opti c al. T his is 
in good agreement with the general finding of Bastian] J2008h that 
r = 0.08 ± 0.03 in different galactic environments. 

The theoretical age distributions all show an initially flat part 
and then a rapid decline. For the lowest mass cut (Mi) the data 
do not allow us to verify this, due to the young age at which the 
Mum(r) line cuts the sample (ti = 10 Myr). The dN/dr of the 
sample with the highest mass cut nicely shows this flat part and the 
decline. The solutions for the mass limited age distributions can be 
approximated within ~10% accuracy by 

* A {rT^H^) ( - 2T/C) (30) 

where t" is the total disruption time of a cluster with an initial 
mass M* and tf is the same for a cluster with an initial mass Mj . 
This approximation only holds for an initial power law index of —2 
at low masses. The approximation in equation ( 1301 ) has the t~ 1//7 
scaling predicted by BL03 for a continuous power law mass func- 
tion, but it falls off faster (exponentially) at old ages, due to the 
truncation in the CIMF. 

A flat part in dN/dr at young ages is a typical result of the 
MDD model (BL03), although BL03 present the age distributions 
for luminosity limited cluster samples, which decline at young ages 
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Figure 5. Ages and masses of M51 clusters (panel [a]) and the resulting mass functions in different age bins (panels [b]-[d]) with the evolved Schechter 
functions (equation |26) over-plotted. The parameters were determined by a maximum likelihood fit of the model to the data (Section [4. l\ . The empirical mass 
functions are approximated by power laws and the resulting index, — f3a t , is indicated in each panel. The mass limited age distributions for different mass cuts 
are shown in panel (e). The vertical positioning of the age distributions with respect to one another is due the fact that each sample has a different lower mass 
limit (Mj ). All mass functions and age distributions are described by a single value of A = 0.07 MQyr~ x , corresponding to T ~ 0.1, i.e. 10% of the star 
formation in M5 1 occurs in star clusters. 



as due to evolutionary fading of clusters. In the prediction 
for dN/dr in this study the evolutionary fading does not enter be- 
cause for each age distribution a lower mass limit is adopted that is 
above the fading line (Fig. [5}- A flat dN/dr was also found for a 
mass limited sub-s ample of clusters in the SMC up to almost a Gyr 
dGieles et alj|2007r) . Indirect evidence for a flat dN/dr for clusters 
in different galaxies follows from the linear scaling of A/ max (log r) 



with t in the first 100 Myr ( Section 12.1.11 and iGieles & Basti"a5 
l2008h - 

Note that there is no evidence for 'infant mortality' from the 
dN/dr in the first ~ 10 Myr, as was reported by B05 from the same 
data set. B05 used slightly smaller age bins and found evidence for 
a drop of a factor of four between the first age bin (< 10 Myr) 
and the second age bin. They also report an enhancement of clus- 
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ters around 60Myr, which coincides with the moment of the last 
encounter between M51 and the companion galaxy (NGC 5195). 
However, there is an artificial age gap between lOMyr and 20Myr 
in the dN/dr, which is a com mon feature when broad-band pho- 
tometry is used to der ive ages dGieles et alj|2005l ; iLee et alj|2005l ; 
Whit more et alj|2007t) . These three effects average out when using 
slightly larger age bins and the details in the dN/dr reported by 
B05 are not separable anymore. The difference between the first 
and the second bin for the sample with the highest mass cut is 
within 1 a still consistent with a decrease of a factor of two or 
three. However, panel (e) of Fig. [5] shows that a flat dN/dr up to 
~ 100 Myr describes the data very well. 

Since this is a quite different interpretation of the data than 
what was concluded by B05, caused by a difference in how the 
data are binned, it is interesting to have a closer look at the dN/dr 
of massive clusters with a method that does not rely on binning. 
One way of doing that is by creating a cumulative distribution 
function (CDF), which can then be compared to different mod- 
els using a Kolmogorov-Smirnoff (K-S) test. In Fig. [6] the CDF 
of the first 100 Myr of the dN/dr of clusters more massive than 
M 3 ~ 6 x 10 4 M Q (Fig[5j is shown (dots). The best-fitting model 
from Section |4~T1 is shown as a full line. The null hypothesis that 
the data have been drawn from this model cannot be rejected. The 
K-S probability is 10%, i.e. the model is within 2a consistent with 
thes e data. 

I Whitmore et alj {2007) claim that cluster evolution is 'univer- 
sal' in the first 100 Myr and that the fraction of disrupted clusters in 
this period is independent of environment and mass. In their model, 
90% of the clusters gets destroyed each age dex, resulting in a r _1 
scaling of the dN/dr for mass limited samples. The CDF of this 
model is shown as a dashed line in Fig. [6] The K-S probability for 
the Whitmore et al. model is 2 x 10 -8 and the hypothesis that these 
data are drawn from a r _1 age distribution can, therefore, be safely 
reject . Even if the disruption fracti on is lowered to 80%( 70%), re- 
sulting in dN/dr oc t~ (, ' 7 (t~°' 5 ) dWhitmore et alj2007l) . the K-S 
probability is 6 x 10 _5 (10" 3 ). 

The fact that the full optically selected cluster population of 
M51 can be described by a CFR that is roughly 10% of the SFR 
indicates that ;/ all stars form in embedded clusters, 90% of them 
are destroyed when the residual gas of the star formation process is 
removed and that this process lasts only a few Myrs. And it seems 
that 'infant mortality' does not affect the optically detected clus- 
ters and that information on the number of embedded clusters is 
needed to estimate the 'infant mortality rate', as was also d one for 
the clusters in the solar neighbourhood dLada & Ladall2003l) . 



5 CONCLUSIONS AND DISCUSSION 

This study considers the early evolution (first ~Gyr) of the star 
cluster mass function, with particular focus on the mass range 
that is available through observations. A Schechter type function, 
i.e. a power law with an exponential truncation (equation [T), is 
used for the cluster initial mass function (CIMF). The use of this 
function is motivated by observational indications that the high- 
mass end of the CIMF of young extr a-galactic clusters is steeper 
than the 'canonical' —2 power law dGieles et al.l I2006a1 ; iLarsenl 
2008). The exponential truncation provides a good description 




Figure 6. Cumulative distribution function (CDF) of ages (^ 100 Myr) for 
the most massive clusters (M > M 3 ~ 6 x 10 ^ MpQ. The dashed and full 
lines show the CDF for the 90% MID model of lwhitmore et al l J2007l) and 
the one presented in this study, respectively. 



of the high-mass tail of the globular cluster mass function (e.g. 
iMcLaughlin & Pudritzll 19961) . 

Empirical cluster masses are generally found to be higher at 
older (logarithmic) ages. This is an observational bias because of 
two effects: the rapid fading of star clusters with age, making it 
increasingly more difficult to see low-mass clusters with increas- 
ing age; also, masses increase due to a size-of-sample effect, since 
longer time intervals are considered at higher logarithmic ages. If 
the CIMF has indeed an exponential truncation, this means that at 
older (logarithmic) ages the steepening is more noticeable than at 
young ages. 

The shape of the cluster mass function at different ages de- 
pends on the functional form of the CIMF and the way clusters lose 
mass due to disruptive effects. Two competing models exist for the 
evo lution of clusters, a mass-dep endent disruption model (MDD, 
e.g. iBoutloukos & Lamers 2003) and a mas s-independent disrup- 
tion model (MID, e.g. lWhitmore et al.i 2007) and both assume that 
star cluster masses are drawn from a continuous power law distri- 
bution with index —2. 

In this contribution a Schechter function is used for the CIMF, 
with a power law index of —2 at low masses and an exponen- 
tial truncation at Al, (equation []}, which is evolved with mass- 
dependent cluster disruption. In summary, it is found that 

• the exponential truncation of the Schechter function is not nec- 
essarily detectable from a small cluster sample, the kind that is used 
to create the mass function of clusters younger than ~ 10 Myr (here 
referred to as CIMF enl p), because in this short time interval there 
are generally not enough clusters sampled above M* ; 

• the age distribution of mass limited sub-samples is flat during 
the first ~ 100 Myr (depending on the mass cut and the disruption 
time-scale) and drops exponentially at older ages. Through a com- 
parison to ages and masses of clusters in M5 1 , it is shown that this 
holds for a sub-sample of clusters more massive than 6 x 1O 4 M . 
The t -1 age distribution, proposed by IWhitmore et all d2007h as 
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the result of a 'universal' cluster disruption model, is ruled out at 
high significance for these data; 

• the mass for which the logarithmic slope of the evolved 
Schechter function is — 2 increases with age as r ' 6 . This scaling is 
similar to the increase of the limiting cluster mass, Mum(r), due to 
the evolutionary fading of clusters. This means that the mass func- 
tion of clusters above the detection limit is approximately a power 
law with index —2 at all ages. 

The MDD model, based on a continuous (nontruncated) power 
law CIMF, predicts that the mass function at old ages should 
have a logarithmic slope of 7 — 1 ~ —0.35 dLamers et al.l 
l2005ah . while empirically derived mass functions at old ages are 
power laws with indices of approximately —2, or even steeper 
(e.g. lde Griis & Anderslkood) . This weakness of the MDD model 
was at the same time a supporting argument for the MID model. 
However, in this study it shown that a power law mass function 
with index —2 at old ages results naturally when a Schechter type 
CIMF is assumed. This is illustrated in Section[4] where it is shown 
that a Schechter function evolved with mass-dependent disruption 
provides an excellent description of the (mass limited) age distribu- 
tion and the mass function at different ages for star clusters in M51. 
The power law mass function of M5 1 clusters gets steeper with age, 
from an index of —2 at the youngest ages to an index of —2.8 at 
~250 Myr. However, simultaneous fitting of models with different 
CIMF and disruption parameters to the ages and masses shows that 
the disruption time depends on mass as tdis oc jVf 67 ± 06 x n i s 
means that the logarithmic slope of the CIMF that is observable at 
~250Myr is smaller than —2.8 (i.e. steeper), since it has already 
been affected by disruption. 

It would be interesting to apply the evolved Schechter function 
to more cluster populations for which age and mass information is 
available. As demonstrated in Section |4~Tl the fundamental param- 
eters defining a cluster population, namely A/, of the CIMF and 
to and 7 of the disruption law (equation I22t can easily be deter- 
mined from lumino sity limited cluster samples u sing a maximum 
likelihood estimate. Parmentier & de Griis] (2008) apply the MDD 
model to the age and mass distributions of LMC clusters and con- 
clude that it is hard to constrain the disruption time-scale. They find 
that models with a long disruption times (> 1 Gyr) are needed to 
describe the mass function, while sh ort disruption times (< 1 Gyr) 
are preferred for the age distribution. lLarse n (2008) showed that an 
exponential truncation in the CIMF around 2 x 10 5 Mq provides 
a good description of the mass function of LMC clusters. Perhaps 
the disruption time-scale of LMC clusters can be better constrained 
when the evolved Schechter function from Section [3~2l is used. 

Most probably the model presented in this paper will not be 
able to e xplain the age distribution of clusters in the Antennae 
galaxies dFall et al]|2005l) . The r _1 scaling of the age distribution 
and the mass-independent nature of the disruption model invoked 
to explain this is very peculiar and this has not been found in other 
galaxies thus far. Since the galactic environment in the Antennae 
galaxies is quite different than that of a quiescent spiral galaxy, 
it could be that clusters suffer from different (additional) mecha- 
nisms that disrupt clusters of different masses equally fast. Some- 
thing along these lines was recently suggested by iRenaud et al.l 
(2008) who studied Antennae like mergers with iV-body simula- 
tions. They show that due to the interaction of the two spiral discs 
there are strong compressive tides that probably induce star and 
cluster formation. Their models show that if clusters or associa- 



tions indeed form in such compressive tides, they can be held to- 
gether for 10 Myr or longer, thereby postponing the dissolution of 
clusters due to gas removal, which in other galaxies seems to de- 
stroy the majority (~90%) of the embedded clusters within a few 
Myrs after formation. Although this is speculative, it is important 
to keep in mind that the cluster system of the Antennae galaxies has 
very different pro perties as compare d to other clust er syst ems (see 
the discussion in iGieles & B astian 20081). IL"arsenl J2008l) showed 
that the CIMF of (normal quiescent) spiral galaxies can be well de- 
scribed by a Schechter function with M* ~ 2 x 10 5 M . It should 
be possible to describe most of the age and mass distributions in 
(quiescent) galaxies with the model presented in this study. 

To validate the correctness of the model presented in this 
study, it would be convincing to 'detect' the turnover in the mass 
function at ages between ~500Myr and ~lGyr. From Fig. [4] it 
can be seen that the depth of observations available at the moment, 
would have to be increased by roughly one dex in mass, or 2.5 mag- 
nitude, to bring the turnover mass, A/to, above the detection limit. 
The improved sensitivity of the new HST/WFC3 camera will al- 
low us to trace the cluster mass function at old ages to the required 
depth to derive a mass function at intermediate ages down to the 
turnover (Mto — 10 4 Mq) for cluster populations at distances of 
5 - lOMpc. 
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